Smoothing a Rock by Chipping 
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We investigate an idealized model for the size reduction and smoothing of a polygonal rock due 
to repeated chipping at corners. Each chip is sufficiently small so that only a single corner and a 
fraction of its two adjacent sides are cut from the object in a single chipping event. After many 
chips have been cut away, the resulting shape of the rock is generally anisotropic, with facet lengths 
and corner angles distributed over a broad range. Although a well-defined shape is quickly reached 
for each realization, there are large fluctuations between realizations. 

PACS numbers: 02.50.Ey, 46.65.+g, 62.20.Qp, 81.65.Ps 
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I. INTRODUCTION AND MODEL 

This work is inspired by a recent experiment of Durian 
et al. [H in which they were interested in the ultimate 
shape of eroding rocks. To investigate this issue quan- 
titatively, they studied the coUisional erosion of square 
clay particles due to their repeated impact with the walls 
of a horizontally rotating plane enclosure. In each such 
collision, chips break off from the particle so that it grad- 
ually becomes rounder and smoother. One might naively 
expect the asymptotic particle shape to be a circle, as 
protruding corners are more exposed and thus likelier to 
get rounded off by this grinding process. 

-N e=7t/2 





FIG. 2: (Color online) Four typical realizations of the cutting 
model for an initial unit square (dashed) after 100 corners 
(circles) . 
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FIG. 1: Schematic illustration of the cutting model after 1, 
2, 3, and 8 events. A chip breaks from a randomly-selected 
corner with the sum of the new deflection angles, 6\ and 62 
equal to the deflection angle 9 of the previous corner. The 
geometry of a single chipping event is shown at the lower 
right. 

A striking result from this experiment is that the 
asymptotic shape of the eroding rock is not a circle [ij. 
To describe this unexpected shape evolution, Durian et 
al. also introduced a simple "cutting" model in which 
the material exterior to a random chord on the object, 
whose length is proportional to the square root of the 
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remaining area, is removed in each cutting event. This 
step is repeated many times until asymptotic behavior 
is reached. Numerical simulations reproduced various 
aspects of the experimental observations and confirmed 
that the asymptotic shape of the particle is not circular 

Q. 

In this work, we investigate an idealized and analyti- 
cally tractable version of this cutting model. A sequence 
of chipping events in our model is schematically illus- 
trated in Fig. [1] The rock is initially assumed to be 
square. In each chipping event, a piece of the rock is 
broken off at a corner. The deflection angles of the two 
newly-created corners sum to the deflection angle of the 
original corner but are otherwise arbitrary. The sides Zi 
and Z2 of a chip are smaller than the respective sides Li 
and L2 of the corner itself (lower right in Fig.[T]), so that 
only a single corner and a finite fraction of its two ad- 
jacent sides are removed in each chipping event. As the 
rock is chipped away, a non-trivial shape is generated 
that is the focus of our interest. 

A convenient geometrical representation for the evolu- 
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tion of corner deflection angles is to view the initial angle 
as a line segment of length 7r/2 (Fig. [3]). Each chipping 
event then corresponds to picking one segment at ran- 
dom and cutting it into two arbitrary size pieces. This 
connection to binary fragmentation allows us to make use 
of well-known results for this latter problem [2] to help 
understand geometrical features of the object as its size 
is reduced by chipping. 
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FIG. 3; Equivalence between angle evolution in chipping and 
fragmentation of a line segment. 



The constraint that the chip breaks off only a single 
corner and a portion of its two adjacent sides ensures that 
final particle size is strictly positive; that is, the initial 
particle is not eroded to nothingness. This fact can be 
appreciated by examining Fig. [2l By the definition of 
the model, a short segment on each of the sides of the 
initial square must remain part of the perimeter of the 
eroding particle in the long time limit. This fact ensures 
that a particle of non-zero size remains in the long-time 
limit. Each of these four segments may be arbitrarily 
small and two such segments on adjacent sides of the 
initial square may be arbitrarily close to each other. In 
the exceptional case where short segments occur in pairs 
at opposite diagonals of the initial square, the area can 
be arbitrarily small, but with vanishing probability of the 
size being zero. 

In the next section, we investigate the evolution of 
the distribution of facet angles by a master equation ap- 
proach when each chipping event bisects the corner an- 
gle. We show that the resulting angle distribution has a 
Poisson form in the variable Ini, where t is the current 
number of corners. In Sec. Illli we study the evolution of 
the angle distribution, as well as the actual shape of the 
object, for the general situation where a chip divides a 
corner angle 9 into two arbitrary angles 9i and 62, with 
= 61+62- For this process, we find: (i) the asymp- 
totic shape of an eroding particle is not round and (ii) 
after many chipping events, the particle is characterized 
by large sample-to-sample fluctuations. In Sec. llVi we 
investigate some natural extensions to more physical cut- 
ting rules and conclude that our main qualitative results 
are robust with respect to these generalizations. We close 
with a brief discussion in Sec. IVl 



II. ANGLE BISECTION 

A. Master equation solution 

We first treat the special case in which a corner angle 
is always bisected in each chipping event. As a result 
of repeated chipping events a non-trivial distribution of 
corner angles develops. Starting with the four right-angle 
corners of an initial square, after one chipping event, two 
angles of magnitude 7r/4 are created, while three right- 
angle corners remain. After a second event, either two 
more 7r/4 angles are created by chipping a right angle 
(this occurs with probability 3/5), or one 7r/4 angle is 
replaced by two new 7r/8 angles (this occurs with proba- 
bility 2/5). 

To determine the angle distribution, it is convenient 
to introduce the integer variable k = — ln2(26'/7r). The 
initial angle 7r/2 then corresponds to fc = 0, and each 
angle halving corresponds to k increasing by 1. We refer 
to a corner with deflection angle corresponding to fc as a 
fc-corner. Let Ukit) be the average number of fc-corners. 
Starting with a square, the initial condition is nk(t = 
0) = 4(5fc^o and the number of corners at time t is simply 
t -\- A. Then the change in nk{t) after one chipping event 
obeys the master equation 



nk{t + 1) - nk{t) = -j—-^ nk-i{t) - 
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Ukit). (1) 



The first term on the right side accounts for the gain of 
fc-corners due to the chipping of one of the (fc— l)-corners 
at time t. The probability of this event is nk-i/{t + 4), 
and each such event increases the number of fc-corners by 
2. Conversely, the second term accounts for the loss of 
fc-corners when one such corner is chipped. 

The system of master equations is recursive, and they 
can be solved one by one. Since n_i = 0, the average 
number of the 0-corners (right-angle corners) satisfies the 
closed equation 



t + 3 

no{t+l) = j---^no{t). 

Iterating this equation, the solution is 

12 



no{t) = 



t + 3 



The average number of 1-corners satisfies 

X t + 3 2 12 
ni(t + 1) = ni(t) H , 



(2) 



(3) 



(4) 



with solution (subject to the initial condition ni(0) = 0) 
, , 24 v-^ 1 

3<j<t+2 

With the solution for ni, the average number of 2-corners 
satisfies 



, t + 3 , , 2 24 

7l2(< + 1) = n2(t) + 

^ ' t + A ' t + At + 3 ^ 3 
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whose solution is 

ni{t) 



48 
iT~3 



E 



3<i<j<t+2 



1 



(7) 



While these exact expression for Uk become progres- 
sively unwieldy as k increases, the asymptotic behavior 
follows easily by noticing that in the t ^ oo limit the 
master equation ([1]) turns into a differential equation 



dnk rik , 2 



(8) 



These equations can also be solved straightforwardly in 
a sequential manner. Using the fact that tiq ^ 12/t and 
rewriting ([8]) as 



d[tnk) 2{tnk-i) 



dt 

we then obtain the solution 



t 



, , 12 (21nt)'= 



(9) 



(10) 



Thus the logarithm of the angle is a Poisson distribution 
with (fc) = 21nt, corresponding to (6*) cx e~*. 

While the final result for nj. is quite simple, we empha- 
size that Eq. pU|) refers to the number of fc-corners aver- 
aged over all possible realizations of the cutting model. 
However, the actual number of /c-corners in a given re- 
alization, defined as TVfe, may differ substantially for 
Ufe = (iVfc). In the appendix, we investigate some of the 
simplest features of the TV^ that illustrate their strongly 
fiuctuating nature. 
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FIG. 4: (Color online) Normalized distributions of corner an- 
gles, Nk/'^f. Nk, versus k for 5 realizations of 10* corners (a) 
and 10^ corners (b). The heavy dashed curve is the properly 
normalized asymptotic expression from Eq. (|10|) . 
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FIG. 5: (Color online) Average and most probable angles for 
a single realization of the cutting process as a function of the 
number of corners. 



B. Simulation results 

We simulated the probabilistic rules underlying the re- 
cursion formula ^ to obtain the distribution for each 
realization. This numerical approach has the advantages 
of simplicity and efficiency, but with the obvious disad- 
vantage that the actual shape of the particle is not ac- 
cessible by this approach. 

The data show that the distribution of the logarithm 
of the angle (actually k — — ln(20/7r)) is close a Poisson 
form in k, as predicted by Eq. (|10p . To compare the data 
with this analytic expression, however, we need to prop- 
erly normalize the latter. Summing Eq. (|10p over all fc, 
one obtains J^k = 12t for the total number of corners, 
whereas the exact result is t -I- 4. To correct for this dis- 
crepancy, we therefore divide the expression in (jlOp by 12 
to compare with the data in Fig. [4l The data are in rea- 
sonable agreement with the properly normalized analytic 
distribution, but it appears that one would have to sim- 
ulate the chipping process for an astronomical number of 
corners to obtain good agreement between the data and 
the asymptotic expression. 

Another important property of the angle distribution 
is the large difference between the average angle {6) and 



the most probable angle ^mp- Because the distributions 
in Fig. [4] are plotted against — ln6', it is clear that {6) is 
much larger than 9^p , which is located at the peak of the 
distribution as shown in Fig. [5l 

Finally, because the natural variable is In 6 the actual 
angle distribution is very broadly distributed. Conse- 
quently, the asymptotic shape of a particle as a result 
of this cutting process will not be circular. A related 
feature is that simulation results from different realiza- 
tions are visually quite different, as might be anticipated 
by the random multiplicative process that underlies the 
chipping process. We will discuss this feature in more 
detail in the next section. 



III. ARBITRARY CHIPPING ANGLES 

A. The angle distribution 

We now study the general situation where a chipping 
event creates two unequal angles. To determine the re- 
sulting angle distribution, we make use of the geometric 
connection between chipping and binary fragmentation. 
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As illustrated in Fig. [3l a corner angle 9 that becomes 
two corners of angles 9i and (with 9i-\-92 = 0), corre- 
sponds to the length-conserving cutting of a segment of 
length 9 into two pieces of lengths 9i and 6*2 . The angle 
distribution in chipping then corresponds to the length 
distribution in the equivalent fragmentation process. 

The length distribution may be solved u&;ing the tech- 
niques from the theory of fragmentation For con- 
venience, consider the scaled segment length x = 29 /tt. 
Starting with a segment of scaled length a; = 1 , the mas- 
ter equation for the length distribution is 



dc{x, t) 
dt 



-c{x. 



P{y, x-y)dy 



c{y,t)F{x,y - x) dy. 

(11) 



Here c{x, t) is the concentration of fragments of length x 
at time t and F[x,y) is the rate at which a fragment of 
size X -|- ?; is cut two pieces of sizes x and y. The first 
term on the right accounts for the loss of fragments of 
size X due to their fragmentation. The total rate of these 
events is F{y, x — y) dy. The second term on the right 
accounts for the creation of a fragment of size y due to 
the breakup of a larger segment of size y. 

In many fragmentation processes 0, S], the breakup 
rate F{x,y) is a homogeneous function of the form 
F{x,y) — {x + y)^^^- That is, the breakup rate of a 
cluster of size x + y depends only on its size and not on 
the size of the two daughter fragments. To make a direct 
connection with cutting, we require A = so that the 
total breaking rate of a fragment is independent of its 
size. In this case, the master equation becomes 



dc{x, t) 
dt 



'C{x, t) 



2 [\iy,t)^. (12) 
Jx V 



This master equation represents the generalization of ^ 
to continuum angles. 

For the initial condition corresponding to a square, 
c(a;, t = 0) = 4:5{x — 1), the distribution of unsealed frag- 
ment sizes at any later time is given by |2| 



c{9,t) = - 



2t 



ttV ln(7r/26') 



e-^h [^JmTL{■K/2e) 



(13) 



TT 



where Ii is a modified Bessel function of order 1. The 
second term on the right-hand side corresponds to the 
probability that there have been no fragmentation events 
up to time while the first term on the right gives the 
scaling part of the fragment size distribution. 

From the variable combination in this first term, we 
find that the characteristic angle has the time dependence 
9 ~ e~*. Furthermore, from the asymptotic form of the 



Bessel function [4|, the distribution of the logarithmic 
angle has a stretched exponential tail 



c(x, t) 



with — In 6* being the natural variable of the system. As 
in the case of symmetric chipping (angle bisection) , the 
general chipping process leads to a broad distribution of 
angles. This result again suggests that the asymptotic 
shape of the particle is not circular. 



B. Shape Evolution 

Area and Perimeter Distributions. Two fundamental 
characteristics of an object's shape are its area and its 
perimeter. Starting with a unit-area square, the resulting 
area and perimeter distributions become smooth, sharply 
peaked about their average values, and visually indepen- 
dent of the number of corners N when > 20. The 
support of the area distribution is [0, 1], with a peak near 
0.67. Similarly, the support of the perimeter distribution 
is [2-^/2,4] and the peak occurs at approximately 3.3. 
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FIG. 6; The distribution of (a) area and (b) perimeter after 
50 corners for 10^ realizations. 

An amusing unexplained feature is that a careful 
examination of the data reveals that the first deriva- 
tives of both distributions are actually discontinuous — 
the area distribution at area equal to 1/2 and 
perimeter distribution at scaled perimeter Pscaicd = 

{P - Pmin)/ (Pmax - Pmin) alsO at ^scaled = 1/2. 

Asymmetry and Fluctuations. After a particle has ap- 
proximately 50 corners, a given realization is visually 
close to its asymptotic shape. As illustrated in Fig. [21 
large fluctuations between different realizations arise, so 
that the shape of a single realization has little connection 
to the average shape. 

Each individual interface typically consists of a few 
longer facets that are punctuated by regions with many 
short facets, with a consequent large change in the lo- 
cal tangent. To illustrate this punctuated interface, we 
show the facet length distribution for 10^ realizations 
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FIG. 7: (Color online) Probability distribution of facet 
lengths P(L) versus L for 10^ realizations of 10, 40, and 80 
corners. The steeper curve corresponds to larger TV. 



with N = 10,40, and 80 corners (Fig. The spike 
&i L — 1 corresponds to the initial unit-length facets 
that remain unchipped. The tail for L > 1 corresponds 
to an initial cut that is sufficiently close to the main di- 
agonal of the initial square so that the facet length can 
be greater than 1 — in fact, the maximal facet length is 
^/2. This large-L tail is distinct from the rest of the dis- 
tribution when the number of corners is small. As the 
number of corners increases, the number of short facets 
correspondingly increases, and there is a huge buildup of 
the small-length tail of the facet length distribution. 
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8: Asymmetry ratio i{N) = ^ {R\{N)) / ^ {BjL{N)) 



versus 1/N. Each data point is based on 10 
with A = 10, 20, ... , 1280. 



realizations. 



Finally, we study how the asymmetry of the particle 
evolves during the cutting process. The proper measure 
of asymmetry is through the moment of inertia tensor of 
an object. For the cutting model, this leads to a cum- 
bersome calculation when the number of corners is large. 
We therefore adopt a simpler approach that should reveal 
the same type of information as the inertia tensor. After 
an initial square has been reduced to an object with a 
specified number of corners A^, we determine the x- and 
y-coordinates of each corner about the center of the initial 



square and then compute the mean-square displacements 
of the X- and y-coordinates of all the corners 



1 



N 



i=l 



1 ^ 

1=1 



in each realization. For each realization we then define 
the larger and the smaller of these two mean-square dis- 
placements, 

RliN) =max{X'^{N),Y^{N)) 
RliN) = mm{X^ {N),Y^{N)), 

and then average these maximal and minimal mean- 
square radii over many realizations. Finally, we quantify 
the asymmetry by the dimensionless ratio 



^(A^) ^ ^ {RliN)) {RUN)). 

Thus, for example, a rhombus in which the two corners in 
the x-direction are a unit distance from the origin while 
the other two corners are a distance 1-l-e, ^ — (1 + e)- As 
a function of 1/A^, the data for ^(n) are quite linear for 
A^ between 10 and 1280. Extrapolating the data of Fig. [8] 
to 1/N 0, we infer the value of £_{N ^ oo) w 1.548, 
with a subjective uncertainty of 0.001. 



IV. EXTENSIONS OF THE MODEL 

An unrealistic feature of our cutting model is that each 
corner has the same probability of being chipped. As a 
consequence, corners tend to congregate, as seen in Fig. [21 
In the equivalent fragmentation process, equiprobable 
corner chipping corresponds to an overall breakup rate 
for a given segment that is independent of its length. This 
length-independent breakup rate in fragmentation de- 
marcates the boundary between scaling solutions, when 
the breakup rate grows with segment length, and "shat- 
tering" solutions [3, 3 , when the breakup rate decreases 
with segment length. The shattering solution is char- 
acterized by a finite fraction of the system being trans- 
formed into a dust of zero-length particles that contain 
a finite fraction of the initial length. This singularity is 
parallel to the gelation transition in irreversible aggrega- 
tion. 

Thus a natural question is whether different behav- 
ior arises in the physically more realistic situation in 
which larger protrusions are likelier to be chipped. In 
the language of the equivalent fragmentation process, we 
should study break-up rates with a positive homogeneity 
index — namely, larger fragments are more likely to break. 
To study the role of a positive homogeneity index, we 
considered the extreme situation in which only the most 
susceptible corner breaks in a chipping event. We thus 
investigated the following three extremal dynamics rules: 

1. Chip the corner furthest away from the origin. 
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2. Chip one of the corners on the longest facet. 

3. Chip the corner with the largest deflection angle. 

Each of these chipping rules focuses on some aspect 
of the most prominent non-smooth regions of the object. 
Qualitatively, we find that these three rules all lead to a 
non-circular asymptotic shape of the object. The reason 
for this non-circularity ultimately stems from the strong 
role played by the first few chipping events. The size 
of each chip can, in principle, range from zero to its 
maximum attainable size (see Fig. [Ij. If one of these 
early chips is close to its maximal size, this chip leaves 
an imprint on the object that persists in the long-time 
limit. This property is different than that of curvature- 
driven interface evolution, in which the amount by which 
a curved region of the interface moves is strictly propor- 
tional to the local curvature [al- 




FIG. 9: Geometry of a two-corner chipping event. 




FIG. 10: Average area of a unit square after A'^ cuts versus 
on a double logarithmic scale when the probability of 
2-corner chipping p2 — 0.9. 

Another natural concern about the applicability of the 
cutting model is the restriction to breaking only a sin- 
gle corner and portions of its adjacent sides in a sin- 
gle chipping event. Indeed this rule ensures that the fi- 
nal size of the particle remains non zero as mentioned 
in the introduction. To test the robustness of the cut- 
ting model results to the possibility that more than one 
corner can be chipped away, we studied the situation in 
which a chip could encompass two corners, as illustrated 
in Fig. [5) Specifically, we pick a corner at random; with 
probability p2, the chip includes both this corner and its 
nearest neighbor. With this rule, the restriction that a 



small segment of the initial square must remain as part 
of the boundary of the eroding particle no longer applies. 
Thus it is not obvious a priori that the size of the ob- 
ject will remain non-zero in the long-time limit. Never- 
theless, this more generous two-corner chipping rule still 
leads to a non-zero particle size, as long as the proba- 
bility for single-corner chips is non-zero. Initially, the 
area decreases rapidly as the number of cuts increases. 
However, As the number of corners becomes apprecia- 
ble, later cuts remove only a tiny fraction of the particle 
so that the area eventually saturates to a non-zero value 
(Fig.IinD. 



V. DISCUSSION 

We studied the geometric properties of an idealized 
model for the erosion of a two-dimensional rock by re- 
peated chipping of small pieces. A chipping event is de- 
fined by cutting a small piece from the rock in which 
a single corner and part of its two adjacent sides are re- 
moved. In our model, each corner has the same probabil- 
ity of being chipped. A two basic outcome of this cutting 
model is that there is shape asymmetry in the long-time 
limit. Thus the asymptotic outcome after many chip- 
ping events is not a circle, as was initially observed in 
the experiments and the simulations of Durian et al. 
Another important feature is that there are large shape 
fluctuations between realizations so that the outcome of a 
single event is not representative of the average behavior. 

We determined the evolution of the distribution of an- 
gles from its governing master equation. For the case of 
angle bisection in each chipping event, we found a broad 
and asymptotically Poissonian distribution of angles in 
the variable in Ini, where t is the number of chipping 
events. This behavior appears to be idiosyncratic to the 
case of angle bisection. In the more realistic case where a 
chipping event divides an initial angle into two arbitrary 
angles (with conservation of the total angle), the angle 
distribution has a different behavior that is immediately 
obtained by the exact correspondence between the distri- 
bution of angles in the cutting model and the distribution 
of fragments sizes in the binary fragmentation of a line 
segment [1, 

Finally, it is worth mentioning that because of large 
sample-to-sample size and shape fluctuations after a 
given number of chipping events, the cutting model does 
not give a unique limiting shape. This behavior is in con- 
trast to the class of interface models where the asymp- 
totic shape of a single realization of interface converges 
to a unique limiting shape. Two famous such examples 
are the strictly convex interface between the origin and 
{x ^ 0,y ^ 0) @ and the interface that is generated by 
the partition of the integers Q- It may be worthwhile 
to explore variants of the cutting model that lead to a 
unique limiting shape to take advantage of the highly- 
developed analysis methods available for this type of in- 
terface evolution process. It should also be of interest 



7 



to extend our study to the more realistic case of three 
dimensions, where there is also a highly-developed math- 
ematical literature on limiting shapes [1]. 
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APPENDIX: FLUCTUATIONS IN Uk 

In this appendix we investigate the statistical proper- 
ties of the distribution of A:-corners in greater detail. We 
show that the actual value of the number of fc-corners in 
a given realization, Nk, is generally quite different than 
the average number of fc-corners, Uk- This lack of self 
averaging can be seen by studying the probability distri- 
bution of the random quantities Nk{t) across all realiza- 
tion. While the full calculation is a tedious endeavor, it is 
fairly simple to obtain the number of right-angle corners 
Noit). The master equation for A^o is simpler than that 
for all the other Nk with fc > 0, because iVo can never 
increase in a single chipping event. 

Starting with a square that has four right-angle cor- 
ners, the number of such corners is a deterministic quan- 
tity when t = and t = 1, 

No{0)=A, 7Vo(l) = 3, 

while for t > 1 the number of right-angle corners is a 
random quantity. Let 

n,(t) = Prob{7VoW 

and let us compute Ilj{t) for t > 1. 

We begin with Hsit). To have three right-angle cor- 
ners, each chipping event must not act on any of these 
corners. Consequently, the probability to have three 
right-angle corners satisfies the recurrence 



Tl3{t + 1) = 



t + 1 
t + A 



Solving this recurrence subject to the initial condition 
113(1) = 1 we obtain 



n3(0 



24 



{t + l){t + 2){t + 3)' 



(A.2) 



By similar reasoning, the recurrence for the probability 
to have two right-angle corners is 

n2(^ + i) = ^n2(t) + ^n3(0, (a.3) 



whose solution is 



^ 36(t-l) 
' it + l){t + 2)it + 3)- 



(A.4) 



The probability to have a single right-angle corner satis- 
fies the recurrence 

+ 1) = ^ ni(t) + ^ u,it), (A.5) 



from which 



ni(t) 



12{t- 2) 
(t + l)(i + 2)(i + 3)' 



(A.6) 



Finally the probability that there are no right-angle cor- 
ners can be found by solving the appropriate recurrence 
formula 



no(i + i) = no(i) 



1 



t + A 



HiW, (A.7) 



or from the normalization 110 + 111+112+113 = 1. In 
either case, we obtain 



(t-l)(i~2)(i-3) 



(i + l)(i + 2)(t + 3) 
As a useful consistency check one can compute 



(A., 



(A.l) and recover the result uq — given in Eq. ([3]) 
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